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We address the problem of the fermion species doubling 
on the Lattice. Our strategy is to factorize the fermion dou- 
bles from the action. The mass term of the Dirac- Wilson 
action is changed. In this case the extra roots which appear 
in the action of free fermions in the moment representation 
are independent of the mass and can be factorized from the 
fermion propagator. However the gauge couplings suffer from 
the pathological ghost poles which are common to non-local 
actions. This action can be used to find a solution of the 
Ginsparg Wilson relation, which is cured from the non-local 
pathology. Finally we compare this factorized action with so- 
lutions of The Ginsparg Wilson relation. We find that the 
present is equivalent to the Zenkin action, and that is not 
quite as local as the Neuberger action. 
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I. INTRODUCTION 

Fermion doubles appear in the lattice whenever one 
calculates the first derivative with the finite difference 
method and Fourier transforms it. This is an old problem 
, which afflicts lattice simulations of any quantum field 
theory with light Dirac Fermions. Bosons are healthy 
because the second derivative has no doubling problem. 
This problem can be addressed in a simple 1 dimensional 
lattice. In this case, with a lattice spacing a, the first 
derivative of a fermion field ■(/;, is calculated in the first 
approximation of the finite difference method, 
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The Fourier transform of an infinite and periodic lattice 
is the Brillouin zone, a continuous interval ] — 7r,7r], and 
the functions are transformed as. 
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which we substitute in eqs. ([^) and find for the derivative, 
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The continuum limit is reached when a vanishes. In this 
case it is natural to replace n x/a and k — s- ap, and 
we get the correct behavior. 
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However, one should check all the points where the nu- 
merator sin k vanish in the Brillouin zone. In Fig. |^ we 
see in fact that sin k has two similar roots in the Brillouin 
zone, at fc = and at = tt. This is the source of the 
doubler problem. 
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FIG. 1. In the Brillouin zone, where k is dimensionless, we 
show the functions sin k, 2 sin ^ and 2 tan ^ . 

In the literature many methods to avoid the doubling 
of fermion species on the lattice are presented. A sim- 
ple solution of this lattice problem would be convenient 
for the study of chiral symmetry breaking in QCD and 
for the study of neutrinos in the Standard Model. Wil- 
son 0] included in the action extra hopping terms which 
projected out the doublers from the onset. However this 
was not chiral invariant, and Kogut 1^ proposed the stag- 
gered fermions to improve chiral invariance. Nielsen and 
Ninomiya produced a topological no-go theorem con- 
cerning chiral invariant lattice models with local charge 
conservation and a local action. To avoid these restric- 
tions, there are some very interesting solutions in the 
literature which depart from the simple Wilson action. 
Kaplan Q extends the lattice with an extra dimension, 
and places the physical fermions in a domain wall, 't 
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Hooft g recommends to have continuum fermions with 
lattice bosons. Neuberger |^ redefines the matrix which 
is used in the fermion determinant, and finds a topologi- 
cally nontrivial solution to the Ginsparg- Wilson [0 rela- 
tion. This opens a new perspective to the study of chiral 
symmetry [§|J9|, and the study of the chiral anomaly on 
the lattice. The method presented in this paper belongs 
to the class of non-local actions [|o|-^ . 

We now show a 1 dimensional example where the dou- 
bles are factorized from the propagator. We find it more 
intuitive to place the fermion fields at the vertices of the 
lattice zones (dual lattice) because we are tempted to 
replace the periodic function sink by the anti-periodic 
2 sin I which has a single root. To get this new function, 
which is wider in the momentum representation, see Fig. 
0, one has to narrow the spacing 2a used in the denom- 
inator of the derivative. We will later shift our lattice 
back to the conventional one, the use of the dual lattice 
here is not crucial. The first approximation of the finite 
difference method is now. 
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and we obtain the desired function. We will use the no- 
tation ipn for the interpolated spinor at the site n of the 
lattice, which is simply the mean value of the spinors 
which are placed at the vertices of the zone. 
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We find that the constant term in sense that it is not 
the first derivative term, has a root aX k — -k. This is the 
extra root which was present in the first definition (|l|) of 
the derivative. 

With the interpolated fields (||) and we define a 
free Dirac action in 1 dimension. 
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were used. We can check that both the mass term (one 
double root) and the derivative term (two different roots) 



of the Dirac action (^ have two roots. Thus this action 
complies with the Nielsen and Ninomiya topological the- 
orem . We see that there are still 2 roots , ^ tan ^ = m 
and fc = TT in the Lagrangian (|7|), but only the first root 
is dynamical in the sense that it depends on the fermion 
mass. Nevertheless it is possible to define a propagator 
G which only shows the correct pole in the denominator, 
and no double. 
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where the cos in the numerator, which cancels the ex- 
tra root in the denominator, comes from the interpolation 
of the fermion fields in the definition of the propagator. 
The extra root is replaced by a pole. The propagator (||) 
has the right continuum limit, 
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For a finite spacing a, we can see in Fig. |l] that this 
propagator also has the peculiar property 2tan| > fc, 
which increases the effect of the ultraviolet cutoff - in 

a 

fermions. This is in contradistinction with sinfc < k. 

In the next Section II we define the factorization pro- 
cedure in detail, for a Lattice in d dimensions. In section 
III gauge invariance is implemented. In section IV the 
factorization is adapted to the Ginsparg- Wilson relation. 
The conclusion is in Section V. In Appendix A an equiva- 
lent factorization is studied, and in Appendix B we study 
lattices where fermions and antifermions are in separate 
sites. 



II. FACTORIZATION OF THE DOUBLES IN D 
DIMENSIONS 

We now extend the action (0) to the d dimensional 
case, and factorize the doubles from the propagator and 
the fermion determinant. In this section we study free 
fermions which can be exactly solved in the momen- 
tum representation. However the factorizing procedure 
must be applicable in position space, in order to include 
eventually the gauge fields, and to compute the propaga- 
tor and the determinant with the standard lattice tech- 
niques. At this point we have to leave the dual lattice 
which was used just for the introduction. In order to 
write the interpolated spinors, it is convenient to define 
the class of vectors pi, with components /i; = -1-1 , —1. In 
d dimensions, the interpolated spinors are now extended 
to. 
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To translate the fields back to the standard positions ft 
we apply a shift of —fl'/2 on the action This shift 
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leaves the momentum representation of the action invari- 
ant. The action is now, 

n,p,p' j 

where we use the convention of hermitean Dirac matrices 
7j. The action ( p^ can be written in a matricial form, 

where the sites are labeled by the n, n' , which are d- 
dimensional vectors with integer components. Like the 
Wilson fcrmion action , this action includes extra hop- 
ping terms. In the present case the chiral non-invariant 
terms are all included in the mass term. In eq. ( p^ ) 
the matrices / and Dj are translational invariant, in the 
sense that the matrix elements only depend on the dis- 
tance to the diagonal. The Fourier transform of these 
matrices depends on a single momentum. 
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and from eq. (0) we get the Fourier transforms. 



dj{k) = i^tan^/(fc) , I(k) = J|cos2 ^ . (15) 
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We define the quark propagator G with the Green equa- 
tion. 
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This Green equation ( |16D is an extended version of the 
normal one. Here we replace the usual identity matrix in 
the right hand side by the matrix /. The action with no 
doubles is finally defined. 



where the matrix Af, in the limit of massless quarks is 
chiral invariant. Like the Neuberger matrix |6[|, the final 
matrix D is redefined from the initial action (|12|) . Using 
the delta relations ^ we find the momentum represen- 
tation of the action, 
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We see that the new action has no extra roots and has 
the correct continuum limit. This result can be com- 
pared with the literature. Except for the mass term, we 
recover the Zenkin action |l^ ], which was derived in the 
formalism of a path integral with the Weyl quantization. 



III. COUPLING TO THE GAUGE FIELDS 

We now include the coupling of the fermions to the 
gauge fields, using the Wilson Q prescription. In order 
to have a gauge invariant action, the fermion operators 
and ■0 are usually linked by gauge transformation paths 
when they are in different sites. In the case of a elemen- 
tary link J , 



^ - = g«-B("+j/2).J 
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where Af, are the generators of the gauge group, and the 
gauge-invariant version of ipni^ft^j 



(20) 



In the present case we have to connect ipfi and ip^^ ^ p-pi 

of action (|l^) by a gauge link. Because each of the d 
components of [jl — jl')/2 can be equal to ±1 or 0, this 
gauge link can be constituted by a product 11 17 of up to 
d concatenated elementary orthogonal elementary links. 
Then there are different paths which connect V'a and V's' , 
and it is convenient to use the average of all the possible 
paths, or to choose a single path in a random way. With 
the links included, the / and d terms in eq. (|l^) are 
gauge invariant an so is the new action (pT[). 

Including the gauge coupling, the fermion part of the 
action ( p^ is extended to. 
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In the continuum limit B tends to a vanishing constant. 
This action ( ^l|) is the one that we propose for the study 
of light fermions. 

In the limit of a constant B and using the d relations 
(^ the action ( pl| ) can be diagonalized in the momentum 
representation. 



d'^k 



i>{k)D{k + B)iP{k) 



(22) 



which is the minimal coupling of the gauge fields to the 
fermions. 

Computing either the propagator or the fermion de- 
terminant, is crucial to directly measure chiral symme- 
try breaking. In the lattice, the Green equation for the 
propagator. 
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= /„-,r(nC/)s,«' . (23) 

is perfectly defined and can be solved numerically with 
the standard lattice methods Although / is not 

the identity matrix, its nonvanishing matrix elements 
are very close to the diagonal. The fermion determinant 
\DU\ is also calculable. As for the pure gauge term, the 
Wilson action is unchanged in this case. 

However the pole in the action, which was introduced 
to solve the doubler problem, raises a new difficulty |T^] 
when the gauge coupling is included. This can be under- 
stood when we extract pertubatively the coupling of the 
fermion to a gauge field from eq. (p^). The same result 
can be extracted from the vector Ward identity, and we 
get, 

d 



1 + tan 



2 "-J 



(24) 



Now the pole in the fermion Lagrangian (|T^) will intro- 
duce a pole in the coupling say of a photon to an elec- 
tron. Hence the doublers, which were expelled from the 
free theory, reappear through the gauge independent 
interaction. This problem is common to the non-local 
actions, except for the Bictenholz-Wiese perfect action 



IV. A SOLUTION TO THE GINSPARG- WILSON 
RELATION 

We can address at the same token |M] the coupling 
problem and the Ginsparg-Wilson0^ relation, 

D^5+j5D^2rD-f^D, (25) 

where r is a constant of the order of the lattice spacing a. 
In the case where the action complies with the Ginsparg- 
Wilson relation, Liischer |^ showed that the lattice has 
redefined a chiral symmetry, which includes the chiral 
anomaly. Moreover, in the continuum limit of a ^ 0, 
the chiral symmetry and chiral anomaly of conventional 
field theory could be recovered Here we will follow 
the Chiu-Zenkin prescription to solve the Ginsparg- 
Wilson relation, which has a linear version for the inverse 
Dirac matrix, 



(26) 



A solution of this equation is the sum of an homogeneous 
solution plus the particular one. Any chiral invariant 
matrix is a solution of the homogeneous equation. 
We will use the solution of eq. (|2^) with a vanishing 
fermion mass m. A simple particular solution is — r. 
Our solution is then. 



D = 



1 



1 



(27) 



In this matrix Z?, the poles are factorized, this is a simple 
extension of eq. (|2^) . If we turn off the gauge fields, the 
Fourier transform can be computed. 



D = 



^1 Ej- 13 tan -f 
1 + ^ *!E,7jtan^ 



V. il^j-tanf +r^tan2^ 



1 + ^ V ■ tan 



2 fej 



(28) 



which now has a single root, with no other doubles, and 
no pole in euclidian space. 





FIG. 2. On a 24x24 Brillouin lattice with parameters 
a = r = 1 we respectively show the Dirac scalar component 
(a), and the 172 component (b) of the Fourier transform of 
the free action. 
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It is simpler when r = a/2 and in the case where all 
components kj vanish except for fc^. In this case eq. ( p^ ) 
is equal to (ji7i sin fc^ + i sin^ ^) , which matches the 
Wilson action. See Fig. ^ for a plot of the free action in 
momentum space. 

In the case where the fermion has a mass m from the 
onset, a natural choice for the action consists in using 
the Dirac matrix D of eq. (j2^) where we replace the 
homogeneous solution by the full solution of eq. (23), 
including the mass term. The action is then. 



\ 



J n,n' 



(32) 



in a logarithmic plot as a function oi\fi~ n'\, see Fig. |^. 
We find that in the limit of large distances is approx- 



imately equal to exp{—2\fi 



7/10.3 



), so the interactions 



decrease faster than a rational function. The finite range 
interactions, like the Wilson action, are usually called ul- 
tralocal , but there are no solutions of the Gisnparg- 
Wilson relation with an ultralocal rD . 



D 



I + r J2j lj(^3 + 



(29) 



and in the free fermion case the Fourier transform is, 



D = 



«f J2j I3 t^^^ ^ + "^-^ tan^ ^ + m + rm 
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V. COMPARING WITH THE NEUBERGER 
SOLUTION 

We now return to the relevant chiral limit of m = 0, 
and compare the factorized action with the literature. 
The present action is formally very simple, and it is 
equivalent to the Zenkin action. It is interesting to com- 
pare it with the Neuberger action which is presently the 
most promising solution of the Ginsparg- Wilson relation. 

The solution ( p7|) has two interesting regions in the 
Brillouin zone. In the center we have the single root. 
There the chiral breaking term vanishes and D i ^. 
The other interesting region is the whole boundary of the 
Brillouin zone, where the chiral invariant term vanishes 
and D ~ l/r. Thus we have a whole surface where 1 — 
rD = and this implies that the Chiu-Zenkin condition 
det{l — rD) = is verified. This is a necessary 
condition |16| for the topological M index 



trj^D =index D 



n+ 



(31) 



to be nonvanishing. The Neuberger Q action D^ has 
the same single root, and has 2^* — 1 points ^ (at the 
traditional position of the doublers) where Djv = l/r. 

It is also straightforward to study the locality of the 
action D. In the literature we can find different defini- 
tions of a local action. Nielsen and Ninomiya called local 
an action with continuous first derivative in momentum 
space. The present action D falls in this class. Liischer 
is more restrictive [|l^, and calls local an action of class 
C°° , in this class the interactions in position decrease at 
least exponentially with distance. The Neuberger action 
belongs to this class. The present action does not, except 
for dimension d — 1, because at the usual position of the 
doubles its second derivatives do not exist. In d = 2 di- 
mensions we Fourier transform back to the position space 
the free action D{k), and draw 



ln\D 




100 120 140 
1^ ^'1 



FIG. 3. On a 500x500 lattice with parameters a = r — 1 
and in the case of free fermions we show ln{\D\fi^fii) as a 
function of n — n'. 

The index of our D has been recently computed by 
Chiu who found that D is topologically trivial. 

The action of Neuberger Dj^ complies with the index re- 
lation. Thus it seems that the action Dj\f is the best can- 
didate to simulate continuum QCD, including the axial 
anomaly. Because D has apparently no chiral anomaly, 
it would be interesting to apply it to the study of chiral 
symmetry breaking on the lattice with a single flavor. 

In continuum QCD the Attiyah-Singer index theorem 
is verified, but in the lattice this is still an open problem. 
According to Niedermayer |21|, this is scale dependent. 
If the topological object is larger than the scale of the 
action then the index theorem should be complied with, 
however if the topological object has a smaller scale, the 
index relation is not verified. In Fig. |^ we sec that a 
gauge configuration with finite topological should extend 
itself for more than 20 or 40 lattice spacings, in order 
to be larger than the scale of the interaction. In the 
case of Dj^, the scale of the fermionic matrix is just a 
few lattice spacings, an it was checked that Dj^ complies 
with the index relation even in a small 6x6 lattice pO[ | . So 
it would be interesting to use very large gauge topological 
objects to check if they would produce a finite index for 
the fermionic matrix D. 
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VI. CONCLUSION 

With the aim to simulate euclidian QCD, we find a lat- 
tice action where the extra roots corresponding to dou- 
bled fermions are factorized from the fermion propaga- 
tor and from the fermion determinant. Like the Wil- 
son fermion action [Q , extra hopping terms are included 
in a initial action. The gauge invariant action with a 
conventional chiral symmetry has undesired poles in the 
coupling of a gauge boson to the fermion, but the ac- 
tion with the Ginsparg-Wilson-Liischer [QjH symmetry 
does not suffer from the same pathology. The factor- 
ized fermionic action D turns out to be equivalent to 
the Zenkin ||l^,|l6| action. This action complies with the 
Ginsparg- Wilson Q relation and with the Chiu-Zenkin 
p6| criterium. However it is only weakly local and it has 
been checked that it is topologically trivial, at least in 
small lattices. The study of this action might help to un- 
derstand the index relation on the lattice. We conclude 
that the interest in the present action is just theoretical 
because the Neuberger action [|j is a better candidate to 
the simulation of QCD on the lattice. 
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simple, however the gauge lattice is 2 bigger than the 
lattice of section III. 

We will study in more detail another example, where 
the anti fermion fields ipn are placed in the normal lattice 
and the fermion fields ipn+p/2 ^-re placed in the dual lat- 
tice (or vice versa). For d > 1, the lattice is composed by 
plaquettes with 4 elementary links and the Wilson pure 
gauge action can be used here. However the number of 
links per zone, the number of links which meet in a vertex 
and the link length differ from the case of Section III. In 
the case of Section HI these numbers were respectively d, 
2d, and a , while in the present case these numbers are 
respectively 2'*, 2'^ and aVd/2. In the case of d = 4 the 
number of links is now increased by a factor of 4 and the 
number of links which meet at a vertex is now doubled. 
Applying the factorization procedure to the free fermion 
action, 

^ = (i) ^^^O'^s+f ^^^^ 

n.p j 

d, (fc) = tan ^7(fc) , I{k) ^Ucos^, 

I 

we find the same momentum representation (|2^) in the 
final action. 
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APPENDIX A: OTHER POSSIBLE SOLUTIONS 
WITH DIFFERENT LATTICES 

The factorization of doubles can also be applied to dif- 
ferent lattices, where the fermion and antifermion fields 
are placed in different sites. In these cases, unlike the 
solution in Section II, interpolations are unavoidable in 
the action, and we no longer have to comply with the 
Nielsen and Ninomiya theorem. It is possible to find 
actions without doubles and without the gauge coupling 
pathology. The gauge links are chosen to directly connect 
the antifermion sites with the nearest neighbor fermion 
sites, it is not necessary to concatenate a few links in 
order to reach a fermion from an antifermion. 

An example of a lattice where the same factorization 
procedure works is a alternate lattice where the ^ and tp 
fields are placed at alternate sites. The link structure is 
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